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, In these notes we review recent progress (and, in Section O we announce a new result) con- 

£NJ ' cerning the statistical properties of the spectrum of Wigner random matrices. 
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The general goal of random matrix theory consists in determining the statistical properties of the 
eigenvalues (and the eigenvectors) of N x N matrices whose entries are random variables with given 
laws in the limit N — > oo. 

In these notes, we restrict our attention to ensembles of hermitian Wigner matrices, whose en- 
tries are, up to the constraints due to hermiticity, independent and identically distributed random 
variables. This assumption guarantees that all eigenvalues are real. Most of the results that I am 
fS| , going to present can also be extended to ensembles of Wigner matrices with different symmetries 
(real symmetric and quaternion hermitian ensembles). 

Wigner matrices have been first introduced by Wigner in 1955 to describe the excitation spectra 
of heavy nuclei. For complex systems like heavy nuclei, it is practically impossible to write down the 
precise Hamilton operator. For this reason, Wigner assumed the matrix elements of the Hamiltonian 
to be random variables, and he tried to establish properties of the spectrum holding for almost every 
realization of the randomness. Remarkably, Wigner's intuition was confirmed by the experimental 
data. More generally, the study of the spectrum of Wigner matrices can be considered as a starting 
point for the analysis of the spectrum of systems with disorder. It is believed, in fact, that the 



spectrum of Wigner matrices shares many similarity with the spectrum of more complicated systems, 
such as random Schrodinger operators in the metallic (delocalization) phase. 

Definition 1.1. An ensemble of hermitian Wigner matrices consists of N x N matrices H = 
(hje)i<j,£<N ■ The entries have the form 

h i i = ~7TT^ x i l + ^ot) forl<j<£<N 
v N 



ha = —= xn for 1 < j < N 

where {xjj,Xj£,yji}, l<j<N,j<£<N, isa collection of N 2 independent real random 
variables. We will assume that the variables {xj£, yj,e}i<j<£<N have a common law, with Exje = 0, 
and = 1/2. We will also assume that the variables {xjj}^ =l are identically distributed, with 
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Exjj = 0, and lEx^- = 1. Throughout these notes, we are going to assume sub- Gaussian decay of the 
entries at infinity, in the sense that there exists v > such that 



< oo, 



< oo 



for some (and thus all) 1 < j < £ < N. 

The assumption of sub-Gaussian decay can be relaxed to subexponential decay, but then some 
of the statements have to be modified slightly. From the definition above, we observe that entries 
of Wigner matrices scale as N~ 1 / 2 with the dimension N. This choice guarantees that, as N — > oo, 
the spectrum of Wigner matrices remains bounded (in fact, as we will see below, with this choice 
the spectrum is contained, almost surely in the limit N — > oo, in the interval [—2, 2]). 

Observe that Definition 11.11 does not specify the precise distribution for the entries of the matrix. 
In fact, one of the main point of the results that will be presented below is exactly the fact that they 
apply to a large class of Wigner matrices, and that they are independent of the particular choice for 
the probability law of the entries. The analysis becomes much easier if one assumes the entries of the 
matrix to be Gaussian random variables. The Gaussian Unitary Ensemble (GUE) is the ensemble 
of hermitian Wigner matrices (as defined in Definition II. ip where all entries are Gaussian. In this 
case, it is possible to write the probability density for the matrix H as 



dP(H) = const 



-^ H2 dH 
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where dH = | J dhi m dh\ m dhjj 

l<i<m<N j=l 



is the product of the Lebesgue measure of the entries of H. What makes GUE particularly simple 
to analyze is its invariance with respect to unitary conjugation. If H is a GUE matrix, and U is an 
arbitrary (fixed) unitary matrix, then also U HU* is a GUE matrix. It is possible to show that GUE 
is the only ensemble of hermitian Wigner matrices which enjoys unitary invariance. Because of the 
unitary invariance, for GUE it is possible to find an explicit expression for the probability density 
function pn o£ the N eigenvalues: 
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The Gaussian factor reflects the Gaussian distribution of the entries. The correlations among the 
eigenvalues are described by the Vandermond determinant squared 
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1.2 



Using this explicit expression for the joint probability density function of the iV eigenvalues of GUE 
matrices, one can compute the local eigenvalue statistics. For fixed k = 1, . . . , N, we define the k 
point correlation function 



pP(m, ■ ■ ■ ,Hk) = J dfj,k+i ■ ■ ■ dfi N p N (m, . . . fi N ) 



It turns out that, as N — > oo, the local correlation functions of GUE converge, after appropriate 
rescaling, to the Wigner-Dyson distribution. For any fixed k £ N, 
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as N — > oo. Here 




Psc (E) = { 27V X "- if 1^1 < 2 (L4) 

otherwise 



is the semicircle law, which describes the density of the eigenvalues around E (we will discuss the 
semicircle law in more details in Section [2]). The convergence in (|1.3p is pointwise in E and holds 
weakly in the variables xi, . . . ,X}-. One refers to the left hand side of (jl.3p as local correlations 

(k) 

because the arguments of p N vary in an interval of size 1/N. As we discussed above, the choice of 
the ./V dependence of the matrix entries implies that all eigenvalues are contained in a finite interval; 
this means that the typical distance between neighboring eigenvalues is of the order TV" 1 . Hence, 
to observe non trivial correlations, the relevant length scale is exactly of the order iV -1 . Eq. (jl.3p 
was first proven for GUE by Dyson in [2] using the explicit expression (jl.ip and the asymptotics of 
Hermite polynomials. Although (jl.ip only holds true for GUE, and there is no explicit expression for 
the joint probability density function of the eigenvalues of any other ensembles of Wigner matrices, 
the convergence (II. 3p is expected to hold true independently of the law of the matrix entries. One 
expects, in other words, the local eigenvalue correlations to be universal; we will present a proof of 
universality in Section [H In fact, universality should hold even more generally; the Wigner-Dyson 
statistics is expected to describe the local eigenvalue correlations in a large class of systems with 
disorder (for example, random Schrodinger operators in the metallic phase). 



2 Density of States and the Semicircle Law 



The first rigorous result in random matrix theory has been obtained in [22j, where Wigner proved the 
convergence of the density of states to the semicircle law for arbitrary ensembles of Wigner matrices 
(hermitian, real symmetric, or quaternion hermitian). For a < b, let N [a; b] denote the number of 
eigenvalues in the interval [a; b]. The density of the eigenvalues, or density of states, in the interval 
[a; b] is defined as 

.A/" [a; b] 
~ N\b-a\ ' 

where the factor of iV in the denominator makes sure that, typically, P[ a -b] is a quantity of order one. 
The density of states P\ a -b] is a random variable whose precise value depend on the realization of the 
randomness. Nevertheless, Wigner proved that, as N — > oo, P[ a -b] approaches a deterministic limit. 
More precisely, he showed that, for any 5 > 0, 



lim P 

iV->oo 



M[a; b] 



Nb-a 



ds p sc (s) 



> 5 







where p sc (s) denotes the semicircle law defined in (jl.4p . After considering the limit N — > oo, we 
can also let r\ = 1 1» — a. | — > 0; in this limit the density of state converges to the semicircle law at 
E = (a + b)/2. Hence, Wigner's result can be formulated as 



lim lim 1 



M[E 



1-E + 2 1 



Nr] 



p sc {E) 



> 5 



It is worth noticing that the semicircle law, and therefore the limiting density of states, is independent 
of the choice of the probability law for the entries of the Wigner matrix. It is also important to observe 
that Wigner's result concerns the density of states on macroscopic intervals, that is intervals which 
typically contain order N eigenvalues (for fixed rj > 0, the interval [E— (rj/2); E+(rj/2)] contains order 
N eigenvalues). What happens if we consider the density of states in smaller intervals, namely in 
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intervals whose size shrinks to zero as N — > oo? These intervals will not contain order N eigenvalues. 
However, as long as the number of eigenvalues is large as N — > oo, we may expect the fluctuations 
of the density of states to be negligible in the limit. This is the content of the next theorem, which 
was proven in [9] (extending previous results from [3 18] ) . 

Theorem 2.1. Consider an ensemble of hermitian Wigner matrices as in Def. Ii.il Let \E\ < 2. 

Then, for any K > 0, 



uniformly in N , for all N large enough. This quantitative bound implies, in particular, that 



This theorem establishes convergence of the density of states to the semicircle law on microscopic 
intervals, that is on intervals containing, typically, a constant (N independent) number of eigenvalues. 
From convergence on the microscopic scale, we also obtain convergence to the semicircle law on 
intermediate scales; for any sequence rj(N) > such that n(N) — > and Nn(N) — > oo as N — > oo 
we have 



Note that, if n(N) < 1/N, the fluctuations of the density of states are certainly important, and one 
cannot expect convergence in probability; in this sense, (12.5P establishes convergence to the semicircle 
law on the optimal scale. 

The proof of Theorem 12.11 relies on two main ingredients, an upper bound on the density of states 
and a fixed point equation for the Stieltjes transform of the empirical eigenvalue distribution. 

Upper bound. Consider a sequence r)(N) > (log N) 2 /N. Then there are constants Kq,C,c such 
that 



for all K > Kq and N large enough. 

A similar statement is valid also for smaller intervals, of size 2/N < n(N) < (log N) 2 /N, but its 
proof is more involved (see Theorem 5.1 in [9]). 

To show (|2.6|) . we observe that (using 77 = tj(N)) 





(2.5) 





(2.6) 
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where {^i a }^ =1 are the eigenvalues of H. Hence 




— — — < — Im Tr = — Im N 

Nr] ~ N H -E-irj N ^ H-E-i-n 




(2.7) 



To bound the diagonal entries of the resolvent, we use (for the case j = 1): 
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(1,1) 
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(2.8) 



H - E-i-n 



hu- E-irj - (a, (B - E - in)- 1 *) 
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where a = (h,2i, ^31, • - • , fyvi) G C^ -1 is the first row of H after removing the (1, l)-entry, and B is 
the {N — 1) x (iV — 1) minor of H obtained by removing the first row and the first column. Using 
the spectral decomposition of the minor B, we find, 



H - E -if] 



XI) 



h 



ii 



E-in- i 

^ u l N l^a=l 



(2.9) 



\ a — E—ir] 



where = N\& ■ u a \ 2 and where {X a }^_^ and {Uo}^ 1 are the eigenvalues and, respectively, the 
eigenvectors of the minor B. We conclude that 



Im 



1 



H-E-vq 



(1,1) < 



1 



< 



1 



H ' LUi N 2^a=l \ a -E-iri N 2^a=l \ a -E-ir) 



Since a is independent of the minor B (and therefore of its eigenvectors), we have 
EC a = NEj^^jUaUWii) = NEj2j^u a (j)u a (i) = 1. 



1,0 



hi 



(2.10) 



Moreover, it turns out that £ Q is well concentrated around its expectation. Therefore, up to a set 
with small probability, we find from (12.10P that 



K 2 



\ a —E—ir\ 



Im4fTr 



(2.11) 



N B-E-ir) 



for a sufficiently large constant K > 0. A more precise analysis shows that the measure of the 
excluded set (on which the bound (|2.1ip may fail) is at most exp(— cy/KNrj); see Lemma 4.7 in [9]. 
Using (|2.7p . but for B instead of H, we obtain 



Im 



H-E-irj 



(1,1) 



< 



K 2 



Nrj 



< 



K 2 



Nrj 



where A/fi[a; 6] denotes the number of eigenvalues of -B in the interval [a; b] and where, in the second 
inequality, we used the fact that the eigenvalues of B are interlaced between the eigenvalues of H. 
Analogously, for any j = 2, . . . , N, we find that 



Im 
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H-E-irj 



K 2 



M[E-%;E+% 



(2.12) 
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up to a set with probability smaller than exp(— c\/KNi]). Since, by assumption, Nrj > (logiV) 2 , the 
total measure of the union of these N "bad" sets is bounded by exp(— (c/2) ^/ KNrj{N)). On the 
complementary set, (|2.12p is correct for all j = 1, . . . , N, and thus, from (|2.7p . we find 



M[E-l-E + l] 
Nrj 



< — 



K 2 



Nri 



M[E-l;E + l] 
Nr] 



<K 



which shows the upper bound (|2.6p . 

Stieltjes transform. We define the Stieltjes transform of the empirical eigenvalue measure by 



m N (z) 



1 



1 



N H-z 



-r— 

N ^ u a -z 
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We are going to compare mi\i(z) with the Stieltjes transform of the semicircle law, given by 



It turns out that Theorem 12.11 follows, if we can show that the difference mjsr(z) — m sc (z) converges 
to zero as N — > oo, for all z = E + ir\ with \E\ < 2 and rj > K/N (see the proof of Corollary 4.2 in 
[7]). It is worth noticing that m sc {z) satisfies the fixed point equation 

m sc {z) + 1 = . (2.13) 

z + m sc (z) 

It turns out that this equation is stable away from the spectral edges E = ±2. Hence to show that 
|mjv(^) — Tn sc (z)\ < CS with high probability, it is enough to prove that mjv(z) is an approximate 
solution of (j2. 13j) . that is that 

1 

m N (z) + 



z + m,N{z) 
with high probability. 

To establish (|2. 14[) . we use again the expression (|2.9p : 



< 5 (2.14) 



1*1 1 N 1 
mN(z) = N^ H^~z^ j) = JV ^ 7 1 x ^ ® ( 2 - 15 ) 

where = JV|a^ • ui j) | 2 , = (fyi, . . . , fyj-i, fyj+i, %iv) G C^ 1 is the j-th row of # without 
diagonal element, and where and are the eigenvalues and the eigenvectors of the minor 
obtained from H by removing the j-th row and the j-th column. From (|2.15p . we obtain 



N 

1 • 

rri]\f(z 

iV 
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N 

1 V - 

JV —z — m N (z) - XU)( Z ) 



1 t 1_A X®(z) 



(2.16) 



Here we defined 



2 + rojv(z) JV ^ (2 + m N {z))(z + mjv(z) + X®{z)) 



X®{z) = - hjj + (^m^_ l( z) - m*(*)) + ^ E ^ 



where 



1 1 1 JV_1 1 

m ^ (*) = ]v^T Tr so) - z = jv^T 2 T(JT 

a=l A a 



is the Stieltjes transform of the minor 

From (|2.16p we obtain (|2.14j) if we can show that X^\z) is small with high probability. The 
first summand in X^\z), —hjj, is of the order JV -1 / 2 and therefore small. The second summand in 
X^\z), the difference {^^-rn'p_ l {z) —rriN(z)), is small (with probability one); this is a consequence 
of the interlacing between the eigenvalues of and H. The main difficulty consists therefore in 
showing the smallness of the random variable 

1 N ~ l t _ 1 
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To show that Y^\z) is small, with high probability, we use the upper bound (|2.6p . which guarantees 
that the denominator in (|2. 17j) can only be large for a small number of a's; details can be found in 
Section 6 of [9]. 



3 Derealization of Eigenvectors of Wigner Matrices 

As a first application of Theorem 12. 11 we show derealization of the eigenvectors of Wigner matrices. 
Given a vector v G with £ 2 norm equal to one, we say that it is completely localized if one of its 
component has size one, and all other vanish, v is called completely delocalized, if all its components 
have the same size (namely TV -1 / 2 ). To distinguish between localized and delocalized vectors, we 
can compute the HP norm of v, for p > 2. If v is completely localized, its £ p norm equals to one, for 
all 2 < p < oo, and for all N. On the other hand, if v is completely delocalized, 

||v|| p = N 2 p 

and therefore converges to zero, as N — > oo (for all p > 2). The next theorem was proved in [9], 
extending results from [7j [8] . 

Theorem 3.1. Consider an ensemble of hermitian Wigner matrices as in Def. Let \E\ < 2, 

K > 0, 2 < p < oo. Then there exists c, C > such that 

3v : ffv = fjcv, \{i-E\< ||v|| 2 = 1, ||v|| p > MN~^ + p^J < Ce~ cVM (3.18) 

for all M > 0, and N large enough. 

If we take M > (logiV) 2 , the statement can be modified as follows. For every tz > 0, and 
2 < p < oo, there exist c, C > such that 

'3v : ffv = fiv, \fi\ < 2 - k, ||v|| 2 = 1, ||v|| p > MN~~i + p) < Ce~ cym (3.19) 



for all M > (logN) 2 , and all N large enough. 

The interpretation of (|3. 18|) and (|3.19p is straightforward; eigenvectors of Wigner matrices are 
completely delocalized. Up to the parameter M > 0, which tunes the probability, all components of 
the eigenvectors have the same size. 

The proof of f)3. 19|) is a simple application of Theorem 12.11 ( (|3.18p requires a little bit more 
work, but the main ideas are the same). For p = oo, we need an upper bound on the components 
of an eigenvector v = (vi, V2, ■ ■ ■ , vn)- To bound, for example, v±, we write v = (vi,w), with 
w = ( v 2, ■ ■ ■ ,vn)- From Hv = pv, we find 

w = —vi(B — /u) _1 a 

where a = (/121, • • • , frjvi) is the first column of H, after removing the first component, and where B is 
the (N — 1) x (N — 1) minor of H obtained by removing the first row and column. The normalization 
condition v \ + w 2 = 1 implies that 

1 1 1 



N l^a=l (A Q -At) 2 N l^a=\ (A a -/ t ) 2 

where, as in Section [2J £ Q = iV|u a • a| 2 , and X a and u a are the eigenvalues and the eigenvectors of 
the minor B. Restricting the sum to those a for which |A a — fi\ < K/N, and using the fact that £ Q 
concentrates around one, we obtain that, with high probability, 

2 < 1 < K 

Vl ~ N \{a : |A Q -fjt\< K/N}\ ~ 
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where, in the last inequality, we use the fact that, by Theorem 12, 1[ the density of states on the 
interval \p — (K/N);/i + (K/N)] converges, with high probability, to p sc (/u) > 0. This shows (|3,19p , 

4 Bulk Universality for Hermitian Wigner Matrices 

Universality for hermitian Wigner matrices refers to the fact that the local eigenvalue statistics 
are independent of the particular choice for the probability law of the entries of the matrix. One 
should distinguish between local correlations close to the edges and in the bulk of the spectrum. 
The correlations at the edges are described by the Tracy- Widom distribution; for GUE this was first 
observed in [21]. In [18], Soshnikov established edge universality of the Tracy- Widom distribution 
(recently, Tao and Vu found a new proof of edge universality in [20]). In these notes, I will restrict 
my attention to universality in the bulk of the spectrum; to show bulk universality, we need to prove 
the convergence of the local eigenvalues correlations to the Wigner-Dyson distribution (j 1 . 3 [) . 

A first partial result towards bulk universality for hermitian Wigner matrices was obtained by 
Johansson, who proved, in [15j . convergence to the Wigner-Dyson statistics for ensembles of Wigner 
matrices with a Guassian component. Johansson considered Wigner matrices of the form 

H = H + t 1/2 V (4.20) 

where Hq is an arbitrary Wigner matrix, V is a GUE matrix, independent of Hq, and t > measures 
the strength of the Gaussian component. Strictly speaking, in order for H to obey Definition ll.il we 
should renormalize it by a t dependent factor to make sure that the variance of the entries remains 
constant. However, this rescaling only affects the eigenvalue statistics in a trivial way, and therefore 
in the following we will work with the matrix H, as defined in (14.201) . 

It turns out that H emerges from Hq by letting every entry evolve through independent Brownian 
motions up to time t. This stochastic evolution of the matrix entries induces a stochastic dynamics 
for the eigenvalues, known as the Dyson Brownian motion. The Dyson Brownian motion can be 
described as a flow for the density of the eigenvalues. The joint probability density function of the 
N eigenvalues of H = Hq + t}l 2 V can be written as 



pnA x ) = J d yqt(x;y)pN,o(y) (4.21) 



where x = (x±, . . . , xjy), y = (yi, . . . , yw), pn,o is the joint probability density function of the 
eigenvalues of the initial matrix Hq, and qt is the kernel describing the Dyson Brownian Motion 
(describing, in other words, the addition of a GUE matrix to Hq). The kernel qt can be computed 
explicitly: 

N N ' 2 A(x) 
{2%t) N l* A(y) 

with the Vandermonde determinant A(x) defined in (|1.2j) . From (|4.2ip . we can express the /c-point 
correlation function oip^t as 



ft ( x;y ) = ^4 det u-Sto-vtffttsN (422) 



where 



pP t (x 1 ,...,x k ) = j dyql k \xi,...,x k ;y)pN,o(y) (4.23) 



q[ k \x 1: . . . ,x k ;y) = I g t (x; y) dx k+1 . . . dx N 

J (4.24) 
(N -kV 

= — — ^{Kt^ixux^y))^^ 



8 



with 



N 

K t ,N(u,v;y) 



(2m) 2 (v - u)t 

N 

w-yj 



f dz [ dw ( e -N(v-u)(w-r)/t _ 1) TT 

h Jr fJi 



w — r + z — u 

N ^ 



t_ Vj r \ c N(w 2 -2vw-z 2 +2uz)/2t 



w-r\ ' N (w - yj)(z - yj) 



Here 7 is the union of two horizontal lines K9s- > s — i5, IR 9 s — > —s + i5, for some 5 > 0, and T 
is a vertical line R 3 s — > k + is and r £ M is arbitrary. 

From (|4.23p . we observe that universality for matrices of the form (|4.20p can be proven by showing 
that 

1 Kt,* ( E + + y) -> - ^ (4.25) 



Np t (E) MV V N Pt (EY Np t (E)' J J ir(x 2 - Xl ) 
for all y in a subset of WL N whose p^o-measure tends to one, as N — > 00. In (I4.25p . 



1 E 2 

Pt{E) 



2vr(l + t) 2 y 4(1 + i) 2 

is the rescaled semicircle law (the rescaling is needed because the variance of the entries of H grows 
with t). 

To show (|4.25|) . Johansson writes 
1 



Kt, N [E,E+ — , -,y 



NpAE) t,JV V ' NpAE)' 

r a r a ( 4 - 26 ) 
= N l^l^l h N (w)g N (z, W )e N ^-M^ 



with 



In(z) = \-{z 2 - 2uz) + — ^log(z - yj ) 



9n{z,w 



2V ' N 

3 

1 L„ „,„„.! 1 \- Vi 



-\w — r + z — u] — —. r ^ 



t(w — r) N{w — r) t-r* (w — Uj)(z — 

h N ( w ) = i( e - T ("-'-)/*«-i) 

and he performs a detailed asymptotic analysis of the integral (|4.26p . He finds the saddle points of 
the exponent /jv(w) — /n(z), he shifts the contours 7, T to go through the saddles, he shows that the 
contributions of the saddles gives exactly the sine-kernel, and that the contributions away from the 
saddles vanish as N — > 00. Following this strategy, Johansson proved convergence to the sine- kernel 
for ensembles of the form (14.20p . for arbitrary fixed t > 0. 

To prove universality for general ensembles of hermitian Wigner matrices, we would like to take 
t = in (|4.20p . As an intermediate step, we may ask what happens if one chooses t = t(N) depending 
on N, so that t(N) — > as N — > 00. The algebraic identities discussed above are still valid, and 
universality follows again by showing that the r.h.s. of (|4,26p converges to the sine-kernel as N — > 00. 
The difference is that now the parameter t entering the definition of /jv(z), gN(z,w), h^{w) depends 
on N and vanishes as N — > 00; this makes the asymptotic analysis more delicate. It turns out, 



9 



however, that, using Theorem 12.11 it is still possible to perform the asymptotic analysis of the 
integral in (|4.26p and to show convergence to the sine-kernel, as long as t(N) > N~ 1+s for some 
5 > O.Theorem l2.1l is really the crucial ingredient of this analysis; the requirement t(N) > N~ 1+s is, 
in fact, a consequence of the fact that we only have convergence to the semicircle law on intervals of 
size K/N or larger. The next theorem was proven in [3]. 

Theorem 4.1. Let Hq be an ensemble of hermitian Wigner matrices as in Def. Choose a 

sequence t(N) such that t(N) -> as N -> oo and t(N) > N~ 1+s for some 5 > 0. Let V be a GUE 
matrix, independent of Hq. Then, the local eigenvalue statistics of H = HQ + t{N) l l 2 V are such that 

1 (k) ( „ %i „ Xk \ f sm(i:(xi — Xj)) \ 

7pf[E+ - ...,£ + * -Kiet //' ^ (4.27) 



f*JE) N \ N Psc {Ey- N Psc {E)) V W^-^)) JiKijKk 

for all \E\ < 2. The convergence here holds after integrating against a bounded and compactly 
supported observable in the variables x\, . . . , x\.. 

One may hope to show that the small Gaussian perturbation cannot change the local eigenvalue 
correlations by comparing the laws of H and Hq. If the entries of Hq have the probability density 
function h, then the probability density function ht of the entries of Hq + t^^V is given by 

d 2 

h t = e tL h, with L = . 

dx z 

At least formally, this implies that \ht — h\ — tLh. Hence, for small t(N) ~ N~ 1+s , the laws of the 
entries of H and Hq are very close. However, to compare the eigenvalue correlations of H and of 
Hq, we would need to compare N 2 entries; therefore, the total distance between the laws of H and 
of Hq is not small in the limit N — > oo. To overcome this problem, we use a time-reversal type idea. 

Suppose that H is an ensemble of hermitian Wigner matrices, whose entries have the probability 
density function h. We can then introduce a new probability density function 

/ y- 2 r 2 f n T n \ 

h=(l-a + t±--.-±l±-)h (4.28) 

and we can define a new Wigner matrix H whose entries have the law h. The probability density 
function of the entries of H-t = H + t l / 2 V is then given by 

h t = e tL h = e tL (l-tL+ t -4 ±-A-\h. 



Universality for the matrix Ht follows from Theorem l4.1l if t = t(N) > N~ 1+s . Moreover, the density 
ht is now much closer to the initial h compared with e tL h. Formally, we find 

fn+l 

Ih-hl^j-^L^h. (4.29) 

Taking t = t(N) ~ N- 1+s , the r.h.s. can be made smaller than any power of JV by choosing n G N 
sufficiently large. Hence, if we choose n large enough, we can compare the correlations of H with 
the ones of H t , and conclude universality for H (from universality for H t ). From (|4.29|) . it is clear 
that, to choose n large, we need to assume sufficient regularity of h; this explains the origin of the 
regularity conditions in the next theorem, which appeared in [4j. 
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Theorem 4.2. Let H be an ensemble of hermitian Wigner matrices, whose entries have the probabil- 
ity density function h(x) = e~ g ( x ' . Fix k > 1 and assume g G C 2 ( n+1 )(M) for some integer n > k/2, 
with 

2(n+l) 

2 <c(i+x 2 r 



3=1 



/or some m G N. T/zen we /iai>e 



^Pn (e+ AT X \^,...,E+ AT X \^)^det ( M*(zi-Zj)) 



p k sc (EY N \ Np. e (EV"' N Psc (E)J V Mzi-x,-)) A< M < fc 

/or a// < 2 (after integrating against a bounded and compactly function 0(xi, . . . ,Xk)). 

Shortly after this result was posted, Tao and Vu obtained, in [19], another proof of universality 
for hermitian ensembles of Wigner matrices. Their proof uses different techniques, but is also based 
on the convergence to the semicircle law on microscopic scales (more precisely, it is based on the 
derealization of the eigenvectors of Wigner matrices, established in Theorem 13 . X [) . The result of 
Tao and Vu requires the third moment of the matrix entries to vanish, but, otherwise, almost no 
regularity. Comparing the two works, we realized that the two approaches could be combined to 
yield an even stronger result. The following theorem was proven in [5]. 

Theorem 4.3. Let H be an ensemble of Wigner matrices. Then, for every k G N, 5 > 0, Eq G 
(—2 + 5; 2 — 5), and for every compactly supported and bounded 0(x\, . . . , Xf~), we have 



rEo+S r 
lim / dE I dx\ . . . dxk 0(xi, . . . ,Xk) 
I->°°Jeo-S J 



' 7# ( E + , • • • ,E + ** ) - det " 



Psc(E) V Np sc {Ey-' N Psc (E)J V (nfa-Xj)) 



(4.30) 



If, moreover, Ex?- = (xij is the real part of the (i,j)-entry of H), we do not need to average over 
E. In other words, we have, for every fixed \E\ < 2, 



lim / dx 1 ...dx k O{x l ,...,x k ) ttj^t pP (e+ Xl ,E+ Xk 
N-+ooJ Psc{E) V Np sc {E) Np sc (E) 

I a a r>f f sm(n(xj - Xj)) \ 
= / dx\ ... dx k U(xi, ...,Xk) det — — ; -— 



(4.31) 



l<i,j<k 



It is interesting to note that the results on universality we just discussed do not generalized easily 
to ensemble of Wigner matrices with different symmetry (real symmetric and quaternion hermitian 
ensembles) . The reason is that the expressions (|4.22p and (|4.24j) are based on a formula due to Harish- 
Chandra to integrate over the unitary group and, unfortunately, there is no analogous formula for 
integrating over the orthogonal (or the quaternion unitary) group. Also the result obtained in [19] 
by Tao and Vu can only be applied to ensembles of real symmetric or quaternion hermitian Wigner 
matrices if the first four moments of the matrix entries match exactly the corresponding Gaussian 
moments. To establish universality for real symmetric or quaternion hermitian ensembles of Wigner 
matrices, we developed new techniques, based on the introduction of a local relaxation flow, which 
approximates Dyson Brownian Motion, but is characterized by a faster relaxation time. In this case, 
we only prove universality after integrating the variable E over an (arbitrarily small) interval (as in 
P~3"0]) ). Details can be found in [10]. 
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5 Extensions and an Application of Universality 



In [T7j, Peche shows universality (pointwise in E) for ensembles of complex sample covariance ma- 
trices extending the approach outlined in Section [H In , we prove universality for ensembles of 
(real or complex) sample covariance matrices (but only after integration over E in a small interval). 

In [T21 H2]) Erdos, Yau and Yin extend the results presented above to so called generalized 
Wigner matrices. Up to the symmetry constraints, the entries of generalized Wigner matrices are 
independent, but not necessarily identically distributed. Instead, one assumes that the variances 
1 < i,j ' < N are such that c\ < No~ij < C2. 



a 



e 4, 



In |14| . the same authors show that the eigenvalues of generalized Wigner matrices are localized 
with very high probability within distances of order (log N) a /N from the position where the semicircle 
law predicts they should be. 

More related results are also presented in the review paper [3]. 

To conclude these notes, we mention a new result concerning the average density of states for 
hermitian Wigner matrices. Theorem 12.11 establishes the convergence of the density of states to the 
semicircle law on the microscopic scale. On shorter scale we cannot expect the density of states to 
converge in probability, because its fluctuations are certainly important. Nevertheless, we may still 
ask whether on these extremely short scales the expectation of the density of states, known as the 
average density of states, converges. The next theorem (see [16]) gives a positive answer to this 
question. 

Theorem 5.1. Let H be an ensemble of hermitian Wigner matrices. Assume that the probability 
density of the matrix entries has the form h{x) = e~ 9 ^ and satisfies the bounds 



Hp) 



< 



l 



(i + c P 2 y/ 2 

Then, for any \E\ < 2, we have 



hg"{p) 



< 



(i + c P 2 yi 2 



for some a > 6. 



• EJ\f 



E-—-E + — 

2N' 2N 



Psc(E) 



(5.32) 



as N —7- oo, uniformly in e > 0. In other words, we have 



lim limsup E 



Af[E-%;E + 
Ne 



^ = lim liminf E l - ^ 

7V->oo e^O Ne 



Psc(E). 



Universality (in the form (|4.3ip ) implies that 



E 



N[ E - m' E + m) 



l(\x\ < k/2) n 1 )/^ x . 

~ Pn( E + Jj) -> Psc(E) . 



(5.33) 



for an arbitrary fixed (N independent) k > 0. Hence, to prove Theorem 15.11 we need to understand 
how (|5.33p can be extended to intervals of size e/N with e going to zero as N — > oo (or even e — > 
before N — > oo). Eq. ()5.32|) follows, if we can show 



- E Im m N (e + i-^j -»• p sc (E) 
uniformly in e > 0. To this end, we establish an upper bound of the form 

Elm m N (E + i— ) £ ( ' ,'Y 

dE V N 



(5.34) 



(5.35) 
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uniformly in e > 0. From (|5.35p . it follows that Elm m^^E + i(e/N)) remains essentially constant 
if E varies within an interval of size n/N, for a small (but fixed) k > 0. This means that 



— E Im m m ( E + i— 
vr V N 

rE4 



~ E- 



N 



7TK 



2N I ( I £ N 

dE Im mAr yE + i— 



E 



E 



-E 

fc L 



E 



arctg 



K K 

2N' + 2N 



arctg 



'N[(JL 



E+ m) 



Psc{E) 



as N — > oo, by (|5.33p . It remains to show the upper bound fl5 .35f> : here we use a Wegner estimate 
from [9] , which implies that the average density of states remains uniformly bounded on arbitrarily 
small intervals. The same techniques used in [9] to show the Wegner estimate can be extended to 
obtain the bound (|5.35|) . Details can be found in [T6] . 
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